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Abstract 

We present a systematic statistical meclianical analysis of the conforma- 



o 

Q \ tional properties of a stiff polyelectrolyte chain with intrachain attractions 

that are due to counterion correlations. We show that the mean- field solution 
corresponds to an Euler-like buckling instability. The effect of the confor- 
mational fluctuations on the buckling instability is investigated, first, qual- 



> 

C^^ \ itatively, within the harmonic ( "semiclassical" ) theory, then, systematically, 

^:^ I within a 1/d-expansion, where d denotes the dimension of embedding space. 

f^ ■ Within the "semiclassical" approximation, we predict that the effect of fluctu- 

^^ ■ ations is to renormalize the effective persistence length to smaller values, but 

not to change the nature of the mean-field (i.e., buckling) behavior. Based on 

J2 \ the 1/d-expansion we are, however, led to conclude that thermal fluctuations 

are responsible for a change of the buckling behavior which is turned into 

'^ ■ a polymer collapse. A phase diagram is constructed in which a sequence of 

^ ■ collapse transitions terminates at a buckling instability that occurs at a place 

Q . that varies with the magnitude of the bare persistence length of the polymer 

chain, as well as with the strength and range of the attractive potential. 
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I. INTRODUCTION 

Recently, the phenomenon of DNA condensation has been subject to renewed interest 
||l||, 0. It appears that the attractive forces generated by the correlated fluctuations of 
or positions of ^ condensed counterions play an essential role in bringing about a sudden 
aggregation of the DNA molecule(s). These interactions are quite intricate since they can be 
described by a pairwise additive potential only at very low concentrations of DNA 0,0. 
For dense DNA mesophases, many-body effects appear to be essential ^ and any theory of 
the DNA condensed state ought to to take them into account. In comparison, the description 
of the solution phase, where pairwise additive interactions make sense, promises to be much 
simpler and more easily within reach. 

Our goal here is to assess some of the fundamental consequences that attractive electro- 
static correlation forces have on the thermodynamic properties of stiff, charged polymers. In 
certain respects this problem has already received some attention in the framework of liquid 
crystalline ordering of polyelectrolye solutions |0]j 0. Here we find it useful to focus entirely 
on single-chain properties. A related approach was introduced recently by Golestanian et 
aii. 

We start with the mean field theory of a stiff polyelectrolyte chain, where thermal fluctua- 
tions are ignored; we find that counterion correlation attractions should induce an Euler-like 
buckling instability of the chain [|10|. Formally, we show that at zero temperature, the buck- 
ling instability is associated with nontrivial solutions of a Shrodinger-like equation which 
is solved for different models of interaction potentials. We determine how Vcr, the critical 
strength of the potential at which buckling sets in, depends on the length of the polyelec- 
trolyte. Below the instability, the chain is straight whereas above the instability point, the 
chain buckles from intrachain attractions. 

In steps, we then treat the effect of the fluctuations on this mean-field picture, first, 
within a harmonic theory that Odijk [|12|, in a related context, has dubbed the "semiclassical" 



picture. This picture is strictly valid only when fluctuation effects are quantitatively small. 



At this level, there are indications that the buckling instability is preserved but that thermal 
fluctuations tend to renormalize the value of the persistence length of the chain, so as to 
make it smaller than the bare value. The harmonic approach breaks down close to the 
buckling instability where conformational fluctuations become large. The predictions can 
therefore only be considered as plausible trends. 

The inclusion of thermal fluctuations of arbitrary strength, however, does not conform 
to the mean-field or the "semiclassical" picture. We propose a new approach that can 
describe properly the thermodynamic properties of the chain for any set of parameters. 



This approach involves a systematic l/d-expansion ||22[ of the partition function, where d 
denotes the dimension of embedding space. We prove that, in general, thermal fluctuations 
completely destroy the (mean-field or "semiclassical") buckling instability and convert it 
into a polymer collapse 0. By varying the bare persistence length, and the strength and 
range of the attractive potential, 

particular, we determine a complete phase diagram. We find that true buckling occurs 
only at zero temperature (or infinite persistence length) and that thermal fluctuations at 
any finite temperature turn buckling into collapse. We do not address the nature of the 
collapsed state (as was done in e.g. ||^); we would have to include consistently into the the 
theory, many-body, non-pairwise additive effects in the electrostatic correlation attractions. 
This, however, is still too difficult. 

As with many concepts in polyelectrolyte theory, the idea that an Euler-like elastic 
buckling instability could lie at the core of the polymer collapse problem has been put forth 



by Manning in relation to DNA condensation IT^]. Manning's calculation, however, does 



not take into account explicitly the attractive part of the interaction potential along the 
polyelectrolyte molecule. The same is true for the recent analysis set forth by Odijk [jl2 



Here we shift the perspective by explicitly considering the contribution of attractive forces 
between polyelectrolyte segments to the onset of a buckling instability. In addition, we show 
how thermal fluctations may transform the buckling behavior into a polymer collapse. 
The exact form of the counterion fluctuation potential in the regime where pairwise 



additivity makes sense, is not of serious concern here. Our theory is vahd for any form. We 
indicate that generally by doing some of the calculations for different model forms of the 
pair potential. However, general arguments based on the analysis of the assymptotic form of 
the counterion fluctuation forces indicate quite strongly 0, |T^ that it is not unreasonable 
to expect that the effective pair interaction potential V{u) is of the generic form V{r) = 
— \vo\('^] e"'"', where r is the monomer-monomer separation, a is roughly the size of an 
effective monomer, vq sets the energy scale for the interaction, and k an inverse screening 
length. We often invoke this form but it is not essential for general conclusions. 

The organization of this paper is as follows. In Section 2 we present the mesoscopic free 
energy of the chain and derive the elastic equilibrium (mean-field) equation for its shape, 
the In Section 3 we investigate numerically the nature of this instability for different model 
interaction potentials. Next we turn to the effect of the fluctuations. The harmonic ( "semi- 
classical" ) theory is considered in Section 4. persistence The analysis based on a systematic 
l/(i-expansion of the partition function, is presented in Section 5. Finally, in Section 6 we 
discuss our results and address the question of possible experimental ramifications of our 
results. 

II. MEAN-FIELD THEORY 

The idea of an interaction-induced buckling instability is really quite intuitive, see Fig. |I|. 
At the point where the elastic energy of bending can no longer compensate for the change of 
the interaction energy due to diminished separations between the interacting segments of the 
polymer chain, the polymer will buckle. Clearly, this buckling depends on the persistence 
length of the polymer as well as on the "strength" of the attractive interactions between its 
segments. 

The starting point for a mean-field description of this buckling instability is the identi- 
fication of a mesoscopic free energy of the self-interacting persistent polymer. We suggest 
here that this free energy may be written as a sum of configurational elastic energy and 



pair-interaction energy for all the monomers: 
Kr r^ ds , r^ fL 






where r(s) denotes a general parametrization of the the polymer chain, s the contour length, 
and |r(s) — r(s')| is the distance in embedding space betweeen two monomers at s and s'. 
Kc is the elastic modulus, which is related to the bare persistence length, ip, via j^ = ip. 
Finally, \^(|r(s) — r(s')|) denotes the monomer- monomer interaction potential assumed to 
be purely attractive (unless indicated otherwise). In what follows we shall remain within 
the framework of the Euler elastic rod theory, ignoring possible non-planar configurations 
pertinent to the Kirchhoffian description |]13| . 



Since we limit ourselves to considering inextensible chains, we ought, in principle, to take 
into account the constraint, dsr{s) ■ dsr{s) = 1, for any point r(s) along the chain. For weak 
undulations, or, in the limit of validity of the mean-field approximation, this constraint can 
be safely ignored [|1^]. In general, however, one has to deal with it appropriately, see Sec. 5. 

In the mean-field approximation we determine a typical configuration of the polymer by 
functionally minimizing the free energy w.r.t. r(s). The result of such a variation is: 






+ jds'¥{\r{s)-r{s')\) = Q, (2) 



where F(|r(s) — r(s')|) = —dV/dr^s) is the local force density. This equation determines 
the typical polymer configuration subject to appropriate boundary conditions. It is easy to 
show that if we consider a polymer with non-interacting monomers subject to the boundary 
conditions that the ends of the polymer do not bend, the solution to the mean-field eqution, 
Eq. (0), corresponds to a straight rod-like configuration. If we imagine that the attractive 
potential between monomers is weak, then deformations away from the rod-like configura- 
tion must be small. Under such circumstances, it is useful to change the parametrization 
of the polymer and to describe the polymer, and solve the mean-field equation, in a coor- 
dinate system which is well suited to describe small deformations away from the straight 



configuration. 

Specifically, for small deviations from a straiglit configuration, tlie polymer cliain can be 
parameterized as r(s) = {z, p{z)), wliere z is tlie direction of tlie axis of tlie molecule and |p| 
is the radial distance from that axis. In this parameterization one has ds = dz\ 1 + ( ^^^ j 

and ^ = — "^^'^ 3/2 ■ With p'{z) = J^ -, the free energy can be written as follows: 

^= f'^dzC[p"{z),p'iz),p{z),z] 
Jo 

= \K,j^dz{p'\z)f{l + {p\z))y' 

+ \l' I' dzdz'^l + (p'(^)) Vl + (P'(^O) V (|r(^) - r(/)|) . (3) 

Jo Jo 



The Euler-Lagrange equation for this free energy reads: 

dC d dC d^ dC 

dp{z) d^dp'iz) ^d^dp"{z) ~ ■ ^^ 

Since here we will be interested only in the limit of a straight rod solution we can linearize 
Eq. (^, deformations, where 

-m - 1< (^-^o^'^' - ^<^'")) - -r-^^^^^^i^ ■ («' 

The first term on the l.h.s. of Eq. (^) stems from the curvature energy. The second is the 
longitudinal stress acting along the deformed rod. The term on the r.h.s. corresponds to 
the transverse bending force |T5|. Both of the last two terms depend on the characteristics 
and the details of the interaction potential. 

By considering the first variation of the free energy at the boundaries, we derive the 
boundary conditions for Eq. (^). Assuming that both ends are free, i.e., that the variations 
of bpi^z = 0, L) and Sp'{z = 0, L) are arbitrary, we find: 

p"{z = 0,L) =0 
lp"'{z = 0, L) - p'{z = 0, L) r dz'V{r{z = 0, L) - r(/)) j = . (6) 

In the case of short range interactions, the r.h.s. term of the Euler-Lagrange equation, 
Eq. d^), can be simplified further. Significant contributions to the integral on the r.h.s. of 



Eq. (^ come only from the points along the polymer for which z and z' are not very much 
apart. Therefore one can develop V(v(^z) — r(z')) in powers oi z — z\ truncating the Taylor 
expansion at the first order term: 

p(z') - p{z) ^ ^(^ - ^') + • • • 
dz 

viz) - viz') = iz- z\ piz) - piz')) ^iz- z') (l, ^) . (7) 

Thus we end up with the following approximation: 

dVi\viz)-viz'))\ _^dVi\u\) P'iz)iz-z') -■^^^\'-''\lp'i,) + ... (8) 



dp{z) du |^_^/|(i + i(p'(^))2 + ...) dz 

To the lowest order in p'iz), the Euler-Lagrange equation reads: 

K,p^'''\z) - 1"^ dz'Vi\z - z'\)p"iz) ^ . (9) 

Eq. (^ is the fundamental mean-field equation describing the shape of the elastic rod in 
the limit of small deformations and short range self-interactions. It is closely related to the 
Schrodinger equation if one introduces the variable u(5;) = p"iz). One finds: 

u"(2) - n^iz)viiz) = with Q^iz) = K-^ f dz'Vi\z - z'\) , (10) 

Jo 

and the boundary condition assumes the form: 

u(^ = 0,L) = 0. (11) 

By definition, Vtiz) is a symmetric function with respect to the midpoint of the rod. By 
integrating uiz) = p"iz), one also gets the part that corresponds to p"iz) = 0, i.e., piz) = 
A + Bz, corresponding to a straight undeformed rod. The field u(z) is therefore a good 
measure of how much the polymer is deformed away from the straight rod configuration. In 
the Eulerian description all u(z)'s are coplanar. 

The mean-field equation, Eq. ([T0|), is not easily solvable except for a limited variety 
of potentials, Viu). Nevertheless, a formal solution can be obtained by introducing the 
following ansatz for uiz), assumed now to lie within a single plane (Eulerian description): 



u{z) =Cri{z) sin {^{z,0)) . (12) 

C is a normalization constant introduced in order to make ri{z) dimensionless. The functions 
ri{z) and $(^, 0) are easily shown to satisfy the following set of equations: 

ri"{z) + n^iz)riiz) - riiz){^f{z, z') = 

{r,\z)^'{z,z'))' = Q. (13) 

The instability point is reached when the following identity is fulfilled: 

1 fL dt 



(14) 



L7o r/2(t) 

which follows from the fact that the non-trivial solution should satisfy the boundary con- 
dition Eq. ( |lT]) which is translated into sin$(L,0) = 0. Eq. (p^) selects the lowest mode 
compatible with this condition. 

Because the function ri{z) still has to be evaluated, this is nothing but a formal statement 
of the properties of the solution of the mean-field equation. The general properties can be 
obtained from a WKB ansatz []T6[. At this level, the comparability r]~'^{z) ~ ^{z), leads to 



a stability limit described by: 

£dzJ\£dz'/3Vi\z-z'\)\ = ./l3IQ7T . (15) 

The necessary condition for the existence of the instability is: 

dz'V{\z-z'\) <0 . (16) 

The properties of this instability (bifurcation) point are the same as in the case of the simpler 
Euler instability where the elastic rod is simply compressed at both ends by a transverse 



force. This latter case has previously been considered by Manning [|TT 



III. MEAN-FIELD SOLUTION FOR DIFFERENT MODEL POTENTIALS 

In order to get a feel for the solutions of the mean-field equation for the shape of the self- 
attracting polymer chain, Eq. ([T0|) , we solve it for three different model interaction potentials, 
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V(r): A finite potential well, an exponential potential, and a counterion correlation potential. 
Approximate analytical solutions can be obtained only for the first two potentials. 
In the case of a finite potential well, 

{—Vo if r < ro 
(17) 
otherwise , 

which leads to. 



n^(z 



Tq < Z < L — Tq 



iz) = K~' f dz'V{\z-z'\) = < 
Jo 



2Voro 

^^^ z < ro (18) 

^^2+-^^^ L-ro<z<L. 



The solutions of the mean-field equation, Eq. (|TU|), with this potential are the angular 
functions and the Airy functions [|T^. Taking into account the boundary condition and the 
continuity of the solution and its derivatives at the discontinuities of the potential, Eq. (ITtI), 
we obtain the following approximate form for the critical magnitude of Vcr'- 

For large enough L, the scaling behavior of Vcr has already been determined by Manning 
pj| . Obviously, the scaling form derived in Ref . [0] is only valid when the condition ^ <^ 1 
is satisfied. 

The next explicitly solvable model is the exponential potential, which can be viewed as 
a generic form of a short range potential. Here 

V{r) = -Voe-''/'^, (20) 

where Vq is a constant. For this potential 

Q\z) = K;' f"^ dz'V{\z -z'\) = ^- ^e-4cosh(^/ro) , (21) 

JO Ac Ac 

where we have displaced the origin of the z axis to the midpoint of the polymer chain, i.e., 
z ^ z + ^. Introducing the variable 2t = z/tq, we obtain the mean- field equation, Eq. (PDj), 
on the form, 

9 



^-(a-2gcosh(2t))y(t) = 0, (22) 

which is equivalent to the modified Mathieu equation with standard parameters a = —^K 

and q = —-^^ exp ^ ||18|. There is only one solution of this equation which satisfies the 

requirements of being both symmetric with respect to the origin of 2;-axis, and finite in the 
limit g — > 0: 

n=+oo , , 
y{t) ~ Y. C2n K2n+s {2^00^^)) , (23) 



where s is the solution of the Hill equation |T^ 



sin^ ^s = -A(0) sinh^ ^^ , (24) 

with Kn{x) denoting the modified Bessel function, and A(0) the Hill determinant for s = 
||T7| . This equation can be solved explicitly only in a limit that would correspond to L/r^ ^ 
1. In this limit A(0) ~ 1 wherefrom: 

sin^ — s ~ — sinh^ — v^? — ^ s ~ i\/a . (25) 

To the lowest order, i.e., for n = 0, the solution of the Schrodinger equation reads: 

y(t)~ir,^(2v/^cosht) . (26) 

The boundary condition at t = ±2^ i\m.?, comes out as 

A-.^(/|(l+e-^)):.0. (27) 

and has to be determined for a = 1 2^ | ■ Only the asymptotic form of the solution can be 
obtained explicitly, as follows: 

not unlike the result for the finite potential well. The details of the interaction potential thus, 
at least in the asymptotic limit, do not appear to matter much. Clearly, the asymptotic form 

10 



derived by Manning |TT| again provides a reasonable description of the point of instability, 
the lowest order deviation from it varying as (^j . 

The last form of the fluctuation potential that we consider is the asymptotic form of the 
effective pairwise additive form of the counterion correlation potential derived in Refs. [H , 



0: 



V{r) = -V,a'[—-^ =-KI^-. (29) 

fluctuation There is no simple analytical result that one can derive for this interaction 
potential even in the asymptotic limit. The numerical solution is, however, revealing enough, 
see Fig. ^ Apparently the finite range of the potential has even less effect on the stability 
limit than in the previous two cases for equal tq. The stability limit is extremely well 
approximated by the WKB ansatz Eq. ([T5| ) which suggests that the critical strength of the 
attractive potential inducing buckling should be inversely porportional to the length of the 
chain squared. 

Comparing the three model potentials one sees that the effect of the finite range of the 
potential is most pronounced for the exponential potential and least pronounced for the 
correlation form, Eq. (p9|). At the mean-field level, the detailed form of the attractive part 
of the intrachain potential is thus of minor importance for a long enough chain. 

IV. FLUCTUATIONS: "SEMICLASSICAL THEORY" 

The mean-field theory analyzed above completely ignores the effects of thermal fluc- 
tuations on the conformational properties of the chain. To include these effects at the 
most primitive level, we now evaluate the partition function of the chain if the mesoscopic 
Hamiltonian is expanded up to second order in the fluctuations around a straight rod-like 
configuration. By analogy with quantum mechanics, Odijk [|I^ has dubbed this type of 



approach the "semiclassical" theory of buckling. The fluctuations are treated very approxi- 
mately in this scheme and as soon as they become large enough, the whole approximation 
breaks down. This happens, of course, right at the instability. Still we will argue that this 

11 



generalization of the mean-field formalism will give us trends concerning the conformational 
properties of the chain close to the buckling transition. 

In order to treat the effect of the fluctuations on the stability properties of the self- 
interacting polymer chain, we first of all write down the free energy of the chain subject to 
small conformational fluctuations. An expansion is performed on the basis of Eq. (^, in the 
limit of p'{z) ^ 1, with the result: 



f-L rL 



J^^IkJ dz{p"{z)f + \l I dzdz' V{\z-z'\) 
Jo Jo Jo 

+ \j^ j^dzdz' [{p\z)f + {p\z')f)v{\z -z'\) + ... (30) 



This harmonic form of the free energy we now write with a new independent variable p'{z) = 
"w{z) on the form, 

J^{w{z),w'iz), z) ~ ^0 + i^c (/'' dz (w'(z))' + I'' dz n\z)w\z)\ , (31) 

with 

J^o = h f I dzdz' V{\z-z'\) , (32) 

JO Jo 



and 



n''{z) =K-^ ("^ dz' V{\z-z'\) . (33) 

Jo 



The propagator for a harmonic action given by the Hamiltonian, Eq. (pi]), can be evaluated 
analytically [|l^] and is determined by quantities characterizing the "classical" solution, as 



calculated via the Euler - Lagrange equation that one can derive from the free energy (or 
"Lagrangian"), Eq. (|3TD. In terms of the local curvature, w'(z), the propagator can be 
derived on the form: 



lC{W{z),z-W{z'),z') = f ■■■ fvw{z) e-/5^(-'(^)'-'(^)^^) 



1/2 (d^/3Fci{w'{z),w'{z'y,z,z')\ 
det ( g^.(,)g^.(,.) lexp(-/3F,,(w(z),w(z);^,z)), 



(34) 
12 



where Fci{W{z)^W{z')]z^z') is the "classical" contribution to the free energy, Eq. (plD, 
evaluated for the w(2;) which is a solution of the Euler - Lagrange equation, Eq. (|10|). The 



propagator can now be written on closed form pO|: 



iri t \ I i\ i\ /sin$(z,2;')\ 

/C(w(.),.;w(.),.)=(-j-^-j-^jx 



/ l3KcCO?,^{z,z' 
exp 



^L {p\z){W{z)f + p\z'){W{z')f) + pK,LW{zW{z') Pl'lff] 



\ 2sm(^{z,z') \r \ J\ \ n r \ j\ \ jj J ^sm^{z,z' 

(35) 

where ^1(2;) = p{z) sin $(2;, z') and V2{z) = p{z) cos $(2;, z') are just two linearly independent 
solutions of Eq. (0). The function p{z) in Eq. (^) is a solution of the Ermakov-Pinney 
equation |pT| , 

p"iz)+n'iz)piz)-^ = 0, (36) 

while $(z, z') can be derived on the form 

All this follows directly from Eq. (|T3|). All other derived quantities can now be obtained 
with the help of this local curvature propagator. 

The density distribution function P[w'(2;)] for the curvature is obtained in a straight 
forward way from the propagator. Assuming free boundary conditions, p"{z = 0,L) = 
w'{z = 0,L) = 0, one obtains the curvature density distribution function as follows: 

P[w'{z)] ~ exp (-^^ (cot $(0, z) + cot $(^, L)) p''{z){w'{z))A . (38) 

The normalization constant is irrelevant as we will be only interested in the average of 
(w'(z)) . Evaluating now the local curvature fluctuations at the midpoint of the polymer 
chain, z = L/2, we are left with 

„ ',rim''\ I<iyy^Ply(^ = L/2)] l (l rt- dt \ 
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where (. . .) indicates thermal averaging. The effect of the fluctuations in the harmonic limit 
can now be assessed as follows. Clearly, at the instability ((w'(L/2)) ) should become large. 
Just how large is difficult to see from Eq. ( pOf ) since the derivation is valid only in the 
limit of small fluctuations. Nevertheless, following Odijk's reasoning [|1^], we claim that the 
instability sets in as soon as the relative fluctuations in the midpoint become larger than 
0(1), i.e., ((w'(L/2))^)p2(L/2) = cq = 0(1). It is difficuh to say more than that without 
actually solving the fundamental equation, Eq. (|36D . Close to the instability point, the 
lowest order contribution to the formal solution of the above equation reads: 



/(L/2){(w'(L/2))=> ~ (,\,L^. ■ («) 

P^cJ^ \2 lJO 2pm), 



Within the WKB approximation this result has a very interesting interpretation. Here the 
above instability limit can be written on the following form, with explicit dependence on 
the interaction potential, 

^ dzJ\ [^ dz'(3V{\z-z'\)\^TxJ]3Kc{l--T^ ^ ^1 . (41) 

If we compare this result with the pure mean-field result, Eq. (|15D, which excludes any effect 
of fluctuations, the instability point is obviously reached when the strength of the potential 
reaches the same value as one finds in the mean-field case, except now the value of the elastic 
constant is renormalized according to: 

^K^ — ^ ^K^ - 2 ^- ^ . (42) 

P^(V2)((w'(L/2))^) ^ ^ 

This renormalization of the elastic constant is obviously fluctuational and is thus linear in 
temperature. Clearly, the above reasoning is not quantitatively valid since we are stretching 
the harmonic theory into a regime where it is not valid. However, one can hope that it bears 
out the correct tendencies for the behavior of this system. 

It thus appears that the "semiclassical " theory of buckling would lead to the same type 
of instability as the mean-field theory, but with the persistence length or, equivalently, the 
elastic modulus taking on a smaller value than the bare value. In other words, thermal 
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fluctuations appear to renormalize the persistence length to a smaller value than the bare 
value. This is exactly the opposite of what happens in the case of purely repulsive segment- 
segment interactions pT| . 

How much of this scenario remains valid in the case where fluctuations can not be dealt 
with as a small perturbation, but are essential to the behavior of the system? The answer to 
this question presupposes the complete solution of the statistical mechanical problem of the 
self-interacting stiff polymer chain. This solution can be obtained only in an approximate 
form such as we derive below. 

V. FLUCTUATIONS: SYSTEMATIC l/L»-EXPANSION 

In order to treat the fluctuations on an appropriate level one has to go beyond the simple 
minded approximate harmonic or "semiclassical" theory we described above. In this section 
we will briefly outline one approach that goes beyond the "semiclassical" theory. We wish, 
in particular, to introduce a program which allows for a reasonably straight forward, ap- 
proximate calculation of the partition function and free energy for a semi-flexible polymer 
whose monomers interact via a pair-potential. The formalism we develop has already been 
applied to describe the conformations and thermal properties of other intrinsically flexible 
materials, including membranes. Thus, the formalism has been used to predict the con- 
formational behavior of fluid membranes |^ and tethered manifolds, with |^|, |2^ and 



without 1^ long-range monomer interactions. In a recent study ^7^ of semi-flexible poly- 
mers with non-interacting monomers {V = 0) an approach, similar to the one introduced 
here, was used. 

For the chains under considerations, we wish to reintroduce the general parametrization, 
r(s), and we wish explicitly to enforce the constraint of "inextensibility" , dsr{s) •9sr(s) = 1. 
The form of the Hamiltonian we shall prefer to use is then, cf. Eq. (^, 
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The partition, itself, is the path-integral over polymer conformations weighted by the Boltz- 
mann weight, exp {—(3T-C), 

Z = /D[r(s)] n S'{dsr{s) ■ dsr{s) - 1) exp (-/3H) , (44) 

where f3 = l/ksT . The functional (5-function guarantees that the integral involves only such 
configurations that satisfy the condition of "inextensibility". 

Two problems complicate the evaluation of the partition function. The first is imposed by 
the functional ^-function and the constraint of "inextensibility" which requires us to include 
in the sum over polymer conformations, r, only those for which the tangent vectors, t = S^r, 
lie on a unit sphere. The second problem which complicates the evaluation of the partition 
function is the fact that, rather than being a simple quadratic (Gaussian) form, the inter- 
monomer interaction potential is, in realistic situations, a complicated, non-local function. 
A systematic way of addressing these problems takes advantage of a "Lagrange multiplier" 
technique. Thus, for instance, one can enforce the constraint of "inextensibility", if one 
introduces an auxiliary field or "Lagrange multiplier", A(s), and adds to the Hamiltonian 
the term, 

1 /■ , .^2 



5ni = -jds\{s){{d,v{s)Y-l). (45) 

Similarly, in order to avoid the complicating non-local term in the pair-potential, 
one can introduce the independent field B = B{s,s'), and make the replacement 
V{{r{s) — r(s'))") -^ V{B\ In order to be able to make this replacement in a system- 
atic way, one must somehow enforce the constraint (r(s) — r(s')) = B. One can do that via 
yet another auxiliary field ("Lagrange multiplier") [Q, |2^ and one is thus led to introduce 
another term in the Hamiltonian, 

m2 = \\ dsd^ g{s,^){{Y[s)-r{^)f -B{s,^)) . (46) 

Given these modifications, the evaluation of the partition function now involves a much 
easier, unconstrained summation over polymer conformations r. The price one has to pay 
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for this simplification is that, in addition to summing over r, one must now sum over X, B, 
and g as weU: 

Z = I D[r(s)]D[A(s)]D[^(s, s)]T>[B{s, s')] exp (-/5(H + 6H, + 8712)) , (47) 

In the expression for the partiotion function, Eq. (^Tj), it is understood that the summation 
over A and g are over contours that begin at —ioo and end at +200. 

It is easy to see that the introduction of "Lagrange multipliers" provides us with an 
expression for the partition function which is quadratic in r and therefore exactly solvable 
as for the integration over polymer conformations. If one fixes \I^ = {A,i?, (7} and expands 
about a particular reference configuration, r = fq, which has the property of minimizing 
7-^2 = 'H + SHi + SH2, i.e., 6H2[r]/Sr\Q = 0, then one finds, after integration, an effective 
Hamiltonian, 

HeK ^ = n2[ro, ^ + ksT^Ti In ( 6{s - s'){ K,dt - d^XOs ) + 2g,{s, s') ) , (48) 

where, 

g^ s') = gis, s') - ^5(3 - s') J ds"igis, s") + gis", s')) , (49) 

and d is the number of components of the vector r or, equivalently, the dimension of em- 
bedding space. If one ignores end effects (by considering a closed polymer, or by enforc- 
ing periodic boundary conditions, say), one can assume that A is a constant, and that 
B{s, s') = B{s — s'), g{s, s') = g{s — s'). It is then possible to perform the diagonalization 
in terms of Fourier modes, so that 

Trln( ... )^ jdsj^H ... ) 

K,dt - Xdl -^ K,q^ + Ag2 (50) 

9c{s,s) ^ g{q) -g{q = 0) , 

where g{q) is the Fourier transform oi g{s—s'). The calculation oi7ie[rQ, \E'] is then, obviously, 
straight forward. 
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What remains in the calculation of the partition function and the free energy, is the more 
difficult integrations over A, B, and g. These integrations can not be performed exactly in 
the general case. If, however, d —>■ cxd, the integrals are completely dominated by the 
contributions from the saddle point, obtained by minimizing w.r.t. A, B, and g. In this 
limit, the exact expression for the free energy of the reference configuration, tq, is therefore 

F[ro] = Fo + (n2K ^] + kBT^Ti In ( K,dt - AS^ + 2g,{s, s') )] , (51) 

where Fq is an unimportant constant. SP implies that the expression is evaluated at the 
saddle point and the Trln-part can be evaluated with the help of Eq. (|50|). For finite d, 
corrections to the saddle point estimate will be of order o{d) and may be calculated via a 
systematic l/d-expansion [^. We shall not do so, being content with the calculation by the 
saddle point method. It is possible to show that this approximation is equivalent to relaxing 
the local constraints, dsr{s) ■ dsr{s) = 1 and (r(s) — r(s')) = B{s,s'), and replacing them 
by the global constraints {dsr{s) ■ dsr{s)) = 1, and ((r(s) — r(s')) ) = B{s, s') [^], ||28 . 



We can now carry out a more quantitative discussion of the properties of semi-flexible 
polymers with pair-wise monomer interactions, which are here taken to be attractive. Within 
the formalism described above, such a discussion can, unfortunately, only be performed for 
simple choices of reference configurations, tq. Here we shall confine our-selves to the choice 



Fq = (se, where e is a one-dimensional unit vector, and C is a "stretching factor" [^, whose 
nature will be described below. 

It turns out that very useful information is contained in the saddle point equations and 
we shall analyze these in some detail. By functionally minimizing w.r.t. A, B{s — s'), and 
g{s — s'), one finds after some manipulations: 

_ 9ro <9ro ,,, rr. f dq f^ , . 

^ " 97 ■ 97 + "^^^^ J 2^K,q^ + \q^ + 2gM) ^ ^ 

gM) = j ds{l~ cos{qs))V\B{s)) (53) 

where V'{z) = dzV . These equations are special cases of more general equations obtained 



by Le Doussal ||2J] and by Palmeri and Guitter ||2^ in their analysis of elastic manifolds 
with long-range monomer interactions. Of these equations, the first, Eq. (0), guaran- 
tees that the constraint, dsr{s) ■ dsr{s) = 1, is satisfied globally, and the third equation, 
Eq. ([5^), takes care of the constraint, (r(s) — r(s')) = B{s,s'). Finally, the second equa- 



tion, Eq. (p5D, determines an effective "self-energy" of the polymer. This "self-energy" may 
be expanded in (even) powers of q, and the expansion coefficients determine contributions to 
the renormalized elastic constants, as may be seen from the expression for the "propagator" 
K~^{q) = l/{Kcq'^ + Xq'^ + 2gc{q)). Roughly speaking, the expansion coefficients tell how the 
non-local interactions modify the parameters involved in a local description of the polymer. 
In particular, V'(B(s)), will, in part, determine a contribution to the total, renormalized 
bending rigidity, as may be seen by analyzing Eq. (|53|) (see further below). 

If, in addition to minimizing w.r.t \1/, one minimizes w.r.t (, so as to determine the best 
choice of configuration in the class of configurations defined by the equation tq = (se, one 



finds, in agreement with Refs. ||2^, [^ 



A(^) = 5jir(g)|,=o = or C = , (55) 



where A'-^^ is a renormalized "Lagrangian multiplier" . li ( ^ 0, the first equation in Eq. 
expresses that straight semi-fiexible polymer is in a stress-free configuration (if the polymer 
were subjected to external stress, the applied stress and A'-'^^ would have to balance). If, on 
the other hand, C = 0^ typical conformations of the polymer will deviate significantly from 
the straight configuration, and we may expect strongly wrinkled configurations to dominate, 
or the polymer to be in a collapsed state (see below). 

It is interesting first to analyze the possibility of having a truly straight configuration, 
with (^ 7^ 0, as the equilibrium configuration of the polymer. Such a straight configuration 
can exist at T = 0, even if the interactions are attractive. However, at non-zero tempera- 
tures we expect that only configurational entropy, which is expected to favor random-walk 
behavior, can prevent the polymer from collapsing and the polymer will, irrespective of 
whether it collapses or not, have equilibrium conformations that deviate significantly from 
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a rod configuration. Indeed, one finds that for C 7^ 0, the integral in Eq. (|52D is disturbed 
by infra-red divergences which can only be removed in the special case Kc = oo. This result 
must be seen as indicating that the straight-rod configuration is unstable. As a consequence, 
in the thermodynamic limit, a phase charaterized by a straight average configuration (the 
ordered phase) exists only in the limit T = (or i^c = oo) and is otherwise destroyed by 
thermal fluctuations. 

Further analysis of the conformational properties of the chain for T > must rely on a 
more detailed analysis of the saddle point equations, Eqs. (^2t)-(^). What can we expect 
from such an analysis ? As we have already indicated, at finite temperatures, we expect that 
polymers whose monomers attract have to compromise between direct interactions which, 
at short scales, favor bending and, at larger scales, favor collapse of the chain in order for 
monomers to be close, and the different effects which counteract these processes, namely 
the initial bending rigidity and conformational entropy. If this picture is correct then we 
must, in agreement with the semiclassical analysis of Sec. 4, expect to find the effective 
bending rigidity, and the persistence length, to decrease significantly, and we must expect 
the theory to signal that collapse of the chain is favorable for strong enough interactions 
between monomers. The quantitative analysis confirms these intuitive considerations. We 
focus on the change of the bending rigidity, which may be obtained by expansion of Eq. (^3|), 

/oo e^ 

ds -rV'{B) . (56) 

-oo 4 

Knowing SK^.^ we can calculate the total, renormalized rigidity as K^^"* = Kc + SK^. 

It is instructive first to consider the situation near T = |^ . At T = we may assume 



^ = 1 and for the interaction potential, V{r) = — -l^exp (— /tr), we find, 

„. = 2^r^e-«-(.. + l) = -i^. (57) 

2 Jo 12 OK 

It is tempting to suggest that this significant reduction in the rigidity could signal that 
the straight rod-like configuration could become unstable even at T = 0. As we saw in 
the previous sections (see in particular Sees. 2-3) such an instability, namely the buckling 
instability, does appear. 
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Now, at finite temperatures, T > (and finite values of the "bare" rigidity, Kc), the 
T = estimate will not describe the conformational properties of the chain very well. We 
must now take ^ = 0, while enforcing the constraint of inextensibility by requiring A^^-* to 
assume a non-zero value. In the following we shall assume that in analyzing gdq)-, it is 
sufficient to retain terms up to and including the fourth order term in the Taylor expansion 
of Qc in powers of q. This is expected to be a valid assumption as long as A and Kc are 
not both small. We then find Qc — (5A)g^ + {6Kc)q'^, where 6X is the contribution to the 
"Lagrange multiplier" from non-local interactions, and A*-^^ = A + (5A. With this assumption 
being made, it is easy to solve Eqs. (E%) and (|5^), with the result: 



BW = ^(?''- + «V"-1)), (59) 



where ^ = vKc /A^^) is a cross-over length which, in the case analyzed here, reduces to 



^ = ^p^Vv^/2- We see that as long as there exists a self-consistent solution for SK^, the 
correlation function B(s — s') = ((r(s) — r(s')) ) behaves as ((r(s) — r(s')) ) oc (^/^p^^)|s — 
s'p oc |s — s'p for small values of s — s', refiecting that the polymer is "straight" on short 
scales. We also see that ((r(s) — r(s')) ) oc (ksT / \^^^)\s — s'\, for large s — s', signaling 
random-walk behavior. 

A self-consistent solution for SK^, may be obtained from Eq. (|56|) after inserting the 
solution to the saddle point equations Eqs. (|58D and (|59|). One derives an integral equa- 
tion whose evaluation, is complicated, for instance, by the cross over between the rigid-rod 
regime and the random-walk regime. In order to obtain the full solution to the problem, 
a numerical study of the integral equation must be carried out. If one is satisfied with 
qualitative/asymptotic results, one can analyze the integral using the method of steepest 
descent. 

The outcome of numerical and qualitative analysis is the following: If one fixes k and 
ip, and determines numerically how i^^^ip (or SKc/Kc) depends on |t>o| one finds results, 
of which the curve shown in Fig. ^: 
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The curve shows that for large values of £^^V^p (small values of IfJ-fCcI) and small values of 
Vq, a solution exists for which i^^^ /ip will decrease {\5Kc{vo)/Kc\ will increase) with |t>o|, as 
it should. For values oi i^^^ip large compared with nip, \{i^^'' —^p)/tp\ = \SKc/Kc\ is found 
to vary linearly with |fo|. This same result is found by a steepest descent analysis, and agrees 
with the T = result displayed in Eq. (0). When \vq\ is increased further, the decrease in 
£^^V^P (increase in \6Kc/Kc\) speeds up. Eventually, for some l^ocl, the rate of change of 
ip^y^p and of \5Kc/Kc\ is "predicted" to become infinitely fast. It is worth noting that |foc| 
is not significantly larger than ksT. Thus, whenever the monomer-monomer interaction 
grows larger than the energy scale set by the thermal energy, the chain will collapse. If one 
then increases Vq beyond Vqc, one finds no solution for 6Kc/Kc which satisfies the demand 
that \6Kc/Kc\ increases as |fo|. This result is in good agreement with the results of a formal 
steepest descent analysis which, for nip -^ 0, predicts that 

6K, = vo^^ X Oil) , (60) 

with a non-trivial relation between 6Kc and fo and the screening length As = 1/n. Physically, 
the latter result is not acceptable and it must be seen as an indication that the theory breaks 
down. In fact, we believe that the solution C = with correlations well described by a random 
walk model will now have to be replaced by a solution characterizing the collapsed state. If 
our analysis is correct, we can therefore conclude that if we investigate the polymer at finite 
temperatures (finite "bare" bending rigidity), we will find that not only will it wish to bend, 
if the interactions between monomers are strong enough, it will also prefer to collapse, in 
order to overcome the entropic penalty. We illustrate these predictions of the conformational 
behavior in the phase diagram. Fig. |^. The interesting feature in this phase diagram is the 
line of collapse transitions which exist for finite values of Kc and terminates at the buckling 
instability point at Kc = c)o. Based on our formal steepest descent calculation, we predict 
that for large enough l/K^, Vc{l/Kc) ~ {l/Kc)~^'^. 
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VI. DISCUSSION 

In the above analysis we explored the connection between buckling of self-interacting 
elastic rods and polymer collapse because of attractive segment-segment interactions. 

This buckling transition is effectively the same as an Euler instability under externally 
imposed compression forces |]T3[. We derived an elastic equilibrium equation whose solu- 



tions determine the state of the elastic rod in the presence of attractive segment-segment 
interactions. In effect, the buckled state corresponds to the bound state solutions of a 
Schrodinger-like equation to which the elastic equilibrium equation is closely related. Ex- 
cept for extremely simple interaction potentials this equation can of course not be solved 
analytically. Nevertheless we find, that the WKB solution quite accurately describes the 
qualitative as well as some of the quantitative aspects of the numerical solution, especially 
with short range potentials. 

Qualitatively, the introduction of thermal fluctuations at the harmonic level, does not 
change the picture of the buckling transition. It nevertheless points to the conclusion that 
conformational fluctuations will renormalize the value of the persistence length. This effect 
is quite well known, if not understood in all its details, in the case of repulsive potentials 
||3 1|| , [ p^ , p3| where the interactions tend to stiffen up the chain. The attractive potential, 
not surprisingly, acts in the reverse direction, thus diminishing the persistence length. The 
harmonic approximation, valid strictly only in the limit of small fluctuations, makes the 
buckling transition, where fluctuations may become prohibitively large, difficult to analyze 
in quantitative terms. We nevertheless argue that it is still there, but displaced towards 
a different point in the parameter space. This displacement is predicted to be linear in 

For unconstrained fluctuations it is difficult to put forth a comprehensive theory. We 
use the systematic l/rf-expansion, which has previously been applied in studies of higher 
dimensional self-interacting manifolds p4| , [^ as a vehicle to build a more general theory of 



the self-attracting polymer chain. On the level of approximation provided by the systematic 
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1/d expansion it appears that buckling in the strict sense of this word is preserved only at 
T ^ or, equivalently, for infinitely stiff chains, K^ -^ oo. At any finite temperature, or 
finite persistence length, the buckling transition is turned into a collapse of the same type as 
already extensively investigated in the case of self- attracting ideally flexible polymers p^ . 

This scenario of course depends on the level of approximation provided by the 1/d- 
expansion in the d ^ oo limit. Usually the variational approach, being non-perturbative, 
does not fare badly; we have confidence that the salient features of the phase diagram for 
the self-interacting stiff polymer chain are not far off from the picture put fourth here. 
The weakest link in our story would appear to be the ansatz tq = (se. It obviously can 
not describe the more realistic toroidal shapes of the e.g. DNA aggregates. But it should 
certainly work fine as long as we are not interested in the detailed structure of the collapsed 
phase but only in the phase boundary. 

At present detailed predictions for experiments are unrealistic. At least the orientation 
dependent part of the interaction should be included in order to describe the nematic nature 
of the condensed state. One thing however we consider to be a robust result of our cal- 
culations: counterion correlation attractions deminish the persistence length. The opposite 
effect of the stiffening of the chain with repulsive intersegment interactions is of course well 
known, though perhaps less well understood |3l||, |^, |^. The effect alluded to here is not 



just the OSF formula with the sign reversed. It has a completely different screening length 
and magnitude dependence than the OSF result. 

The linearized version of this effect is embodied in Eq. p^). Recent experiments on 
stretched DNA in the presence of variable amount of Co(NH3)g"^ [^] clearly show that 



the effective persistence length gets smaller the higher the concentratiuon of the condensing 
agent, that without doubt confers some correlation attraction to the intersegment interaction 
potential. In these experiments the concentration of the condensing agent is too small to 
induce a full blown collapse of DNA, but still, the incipient effects are seen in the smaller 
effective peristence length. Qualitatively this is exactly what one expects from our theory. 
Also the present form of the theory seems to be well suited to describe the effects of the 
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correlation attractions on the elastic extension of the chain. The interplay between collapse 
and stretching seems to be well within the reach of the present formalism and will be pursued 
in all the details in the future p6 . 
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VIII. FIGURE CAPTIONS: 

Figure 1: Schematic representation of buckling. In a buckled chain the effective sepa- 
ration between segments becomes smaller which lowers the free energy of the chain if the 
interactions between segments are attractive. This decrease in the free energy works against 
the increased bending energy thus leading to an instability. 

Figure 2: Critical strength of the interaction potential, in terms of Vcr L^, as a func- 
tion of the length of the rod, L, at different values of the screening parameter tq for the 
three models of attractive potentials: Finite potential well, Eq. ([T7|), exponential poten- 
tial, Eq. (|20|) , and the general correlation potential, Eq. (p9|) . The dependencies have been 
rescaled in such a way that at large L they coincide. 

Figure 3: The solution of Eq. (^61), the relation between the strength of segmental 
attraction and the change in apparent bending modulus, for nip = 0.1. Observe that ^I^V^p 
depends linearly on |fo| for small values of |fo|. There exists a |foc| where there is no longer 
a physically acceptable solution for ^p^V^pd'^ol), implying a loss of stability of the coiled 
configuration of the polyelectrolyte chain. It is worth observing that |foc| is not significantly 
larger than ksT. 

Figure 4: The phase diagram for a semi-fiexible polymer whose bending rigidity is Kc, 
and whose monomers interact via an attractive potential of strength \vo\. Buckling of a rod 
may take place for some |fo| when the bending rigidity is effectively infinite. Collapse of a 
random coil may take place for some |fo(-ft'c)| when Kc is finite. 
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IX. FIGURES: 
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FIGURES 




FIG. 1. Schematic representation of buckling. In a buckled chain the effective separation 
between segments becomes smaller which lowers the free energy of the chain if the interactions 
between segments are attractive. This decrease in the free energy works against the increased 
bending energy thus leading to an instability. 



1.6 
1.4 



1.0 
0.8 



1.2 - 







r„=0.4 




well 


A 








exp 








corr 


\ 
















Unstable region 


- 


■s^ 
















Stable region 








r„=0.2 












Unstable region 


^ 
















Stable region 




2 3 


4 


5 


6 7 


8 9 1 



FIG. 2. Critical strength of the interaction potential, in terms of Vcr LP', as a function of 
the length of the rod, L, at different values of the screening parameter ro for the three models 
of attractive potentials: Finite potential well, Eq. (p^, exponential potential, Eq. ([20|), and the 
general correlation potential, Eq. (|29|). The dependencies have been rescaled in such a way that at 
large L they coincide. 
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FIG. 3. The solution of Eq. (|5q), the relation between the strength of segmental attraction and 
the change in apparent bending modulus, for nip = 0.1. Observe that ip /ip depends linearly on 
\v()\ for small values of \vo\. There exists a |foc| where there is no longer a physically acceptable 
solution for £p Y^p(l^ol)) implying a loss of stability of the coiled configuration of the polyelectrolyte 
chain. It is worth observing that \voc\ is not significantly larger than ksT. 
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FIG. 4. The phase diagram for a semi-flexible polymer whose bending rigidity is Kc, and whose 
monomers interact via an attractive potential of strength |fo|- Buckling of a rod may take place 
for some |fo| when the bending rigidity is effectively infinite. Collapse of a random coil may take 
place for some |wo(-f^c)| when Kc is finite. 
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